Piecewise adiabatic population transfer in a molecule via a wave packet 
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We propose a class of schemes for robust population transfer between quantum states that utilize 
trains of coherent pulses and represent a generalized adiabatic passage via a wave packet. We 
study piecewise Stimulated Raman Adiabatic Passage with pulse-to-pulse amplitude variation, and 
piecewise chirped Raman passage with pulse-to-pulse phase variation, implemented with an optical 
frequency comb. In the context of production of ultracold ground-state molecules, we show that 
with almost no knowledge of the excited potential, robust high-efficiency transfer is possible. 

PACS numbers: 42.50.Hz, 82.50.Nd, 82.53.Kp, 34.50.Cx 



With the advent of technologies that deal with quan- 
tum properties of matter, there is a need for methods of 
controlling quantum systems that are, on one hand, ro- 
bust and conceptually simple, and on the other, flexible 
and applicable far beyond few-level arrangements. To 
address this need, our goal in this work is to utilize the 
intuition and control recipes of simple schemes for ro- 
bust transfer of quantum population based on Adiabatic 
Passage (AP), and apply them to many- level quantum 
systems with complex dynamics. The key components of 
the method are coherent accumulation of amplitudes by 
a train of phase coherent optical pulses, and interference 
of quantum pathways in multi-state population transfer. 

We test our ideas, simulating dynamics in KRb 
molecules, which is of special interest because of the 
prospect for creating deeply bound ultracold polar 
molecules starting from loosely bound Feshbach states 
[TJ [5]. As translationally and vibrationally cold polar 
molecules are a prerequisite for many applications [3], 
their robust high-yield creation is a forefront goal. 

Coherent accumulation has been studied in the context 
of perturbative control [4] and later on in precise spec- 
troscopy using frequency combs [5]. Its use for complete 
population transfers was studied so far twice. One was 
the analysis within the 3 level model of piecewise adia- 
batic passage (PAP) [B]; the other study demonstrated 
population transfer through a wave packet with a coher- 
ent train of pump-dump pulses [7]. Here, we unite the 
two as special cases within one general framework, where 
the combination of adiabatic transfer concepts with a 
coherent train excitation (optical frequency comb) pre- 
serves the robustness of adiabatic transfer, but is appli- 
cable to composite quantum states. 

The concept of piecewise adiabaticity is elucidated 
with two main examples. In piecewise stimulated Raman 
adiabatic passage (STIRAP), robust transfer is achieved 
through a slow variation of the intensity envelope of the 
driving pulse trains. In the other, which we term "piece- 
wise chirped Raman passage" (CRP), robustness is ob- 
tained by slow variation of the excitation phase. We 



demonstrate how both examples can be extended to the 
case of the intermediate state being a wave packet, i.e. 
a multiplicity of intermediate states. In the spectral do- 
main, the transfer is seen as a constructive interference of 
several AP pathways [HE]- Last, we discuss the case of 
transfer without detailed knowledge of the intermediate 
dynamics. We show that by scanning the comb param- 
eters of a train of unshaped pulse-pairs (repetition rate 
and intra-pair time separation) it is possible to achieve 
high efficiency transfer. The intermediate dynamics can 
then be studied using a two-dimensional mapping of the 
transfer efficiency as a function of the comb parameters. 

A traditional AP transfers population from an initial 
to a target state by dressing the quantum system in 
slowly changing fields. If one of the time-dependent field- 
dressed eigenstates coincides with the initial state at the 
beginning of the process, and with the target at the end 

[Till EU IH] , a complete and robust transfer is obtained 
as long as the dressing fields change adiabatically. Ac- 
cordingly, the prescription of PAP [6J can be summarized 
as: For reference, consider any traditional AP, driven by 
slowly varying fields. Then, break the reference AP into a 
set of time intervals r„ , so short that only a small fraction 
of population is transferred between the eigenstates dur- 
ing each interval. Within each interval, the driving fields 
can be replaced by any other field amplitude shape, as 
long as the integral action of each of the fields over r n re- 
mains the same. For example, the smooth reference field 
can be replaced by a train of mutually coherent pulses. 
Last, vary the inter-pulse 'silent' time of the system's free 
evolution. Controlling the silent, free evolution durations 
and the phase of each of the driving fields allows control 
over the arising Ramsey- type interference picture. 

Ref. [5J introduced piecewise versions of STIRAP 
jTTJ [T^] , which transfers population in a three- level sys- 
tem from state |1) into state |3) via the intermediate 
1 2). In the key example, the field was given as a train 
of femtosecond "pump" and "dump" pulse pairs, reso- 
nantly coupling the eigenstates |1) and |2) (pump), and 
the eigenstates \2) and |3) (dump). The pulses - "kicks" 
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FIG. 1: (Color online). Sequences of femtosecond pulses in 
piecewise STIRAP (a) and piecewise CRP (b) between states 
jl) = 35", 1 2} = 4P, and |3> = 5S of Na Q3]. Dashed blue: 
pump field envelope; solid red: dump field envelope; gray: the 
carrier phase in piecewise CRP with ap = old — 0.2. 

- did not overlap in time and the pulse train envelopes 
were slowly varied to achieve PAP, as shown in Fig. [TJa). 

Another conventional AP in a three-level system is 
CRP [TTj . The pump and dump pulses share a similar 
smooth temporal intensity profile; however, both pulses 
are frequency chirped, so that their phases are 4>P,D(t) = 
u P<D t + a { p^ D (t - t ) 2 /2. Like STIRAP, CRP transfers 
population from state |1) into |3) in a robust manner; un- 
like STIRAP, CRP creates transient population in state 
|2). For the piecewise version, all pulses are exactly on 
resonance. The effect of chirp is mimicked by varying the 
carrier phase of each pulse 4>p^ D (n) = 4>p t D( n , i) ~ ^>p nt 
quadratically with the pulse number n: 

4>f D {n) = aPtD {n~n f/2 (1) 

as shown in Fig. |TJb) . Pulse trains with such a piece- 
wise chirp demonstrated piecewise adiabatic following in 
a two-state system [14, . Here we observe that the piece- 
wise CRP in a three-level system is robust with respect 
to the field strengths, the number of pulse pairs, the val- 
ues and signs of ap and otr>, and to additional delays 
between the pump and dump pulse trains. 

In both foregoing examples, the details of the single 
pulse dynamics are not important, only the values of the 
basis state amplitudes after each r„ matter. Based on 
this observation, below we extend the conceptual scheme 
of PAP, generalizing the eigenstate |2) to a multitude of 
states, forming a wave packet. This wave packet under- 
goes complex dynamics between the pairs of pump and 
dump kicks, but if the pulse repetition time coincides 
with a revival, the wave packet (almost) returns to its 
original state by the time the next pair of kicks arrives. 
The kicks, in turn, become rather complex operations 
implemented via an interplay of the shaped femtosecond 
pulses with the intermediate wave packet dynamics. We 
expect that, as long as the action of the kicks on the 
basis states |1), |2), and |3) mimics that of the dressing 
fields in the reference AP coarse-grained over r„, the en- 
forced population transfers in both schemes is insensitive 
to the parameters of the driving fields. Further, high fi- 
delity of the revivals turns out to be unnecessary, as we 
demonstrate later with unshaped pulses. 

Below we study this general idea simulating creation 
of deeply bound X 1 Y, + KRb molecules from a a 3 £ + 
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FIG. 2: (Color online). KRb potentials, basis states, and 
PAP arrangements discussed in the text. 

state. The molecular potentials used in the simulations 
are shown schematically in Fig. [2] Eigenenergies and 
transition strengths are calculated using the algorithm 
FDEXTR [T5] with the data input from Refs. PU E]. 
The excited electronic state is modeled as the LS-coupled 
A 1 !] and 6 3 FI potentials. Although the model is too sim- 
ple to precisely describe molecular dynamics on nano- 
microsecond time scales, it retains the main complexity 
of the problem due to the coupled singlet-triplet evolu- 
tion. We choose |1) to be the vibrational state v = 5 
(E = -157 cm" 1 ) of the a 3 £+ potential (Fig. [2}. To 
trace the fidelity of the method, we also consider the am- 
plitudes of vibrational states around the input state (v = 
0, 1, 12 of the a 3 X + potential). The wave packet |2) is 
composed of up to 20 vibrational states of the LS coupled 
A 1 S-6 3 n potentials with the energies E = 11000 - 11400 
cm" 1 . State |3) is v = 22 (E = -2490 cm" 1 ) of the 
ground X 1 H + potential; the neighboring X 1 Y, + states 
v = 16, ...,28 are included in the simulations as well. 

The system, initially in state |1), is driven by a se- 
ries of mutually coherent pump and dump femtosecond 
pulses. A single pump pulse (110 fs FWHM in inten- 
sity, sin 2 crf-shaped, as shown in Fig. (3^a)), excites an 
A 1 S-6 3 II wave packet, which oscillates on the coupled 
yl 1 S-6 3 II potential surfaces with the singlet and triplet 
vibrational periods of the order of 800 — 900ps. After 
several vibrations the dynamics become irregular. How- 
ever, at certain times, which correspond to two-surface 
revivals, the wave packet rephases with a fidelity reaching 
~ 0.9. In our simulations, high-fidelity revivals were ob- 
served, roughly, every 80 ps; the particular revival time 
is sensitive to the spectroscopic data used as input. Fol- 
lowing its excitation, the wave packet can be dumped 
into the single vibrational level |3). The shape of the 
femtosecond dump pulse is found by requiring that the 
wave packet before dumping overlaps well with the wave 
packet that would have been excited by the time-reversed 
dump, in a manner similar to the prescription of Ref. [7]. 
The thus- found dump field is shown in Fig.J^b). 

The simulation results for both piecewise STIRAP 
and piecewise CRP are presented in Fig. |4j For STI- 
RAP, shown in Fig. |ija,b), the field consists of 200 
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FIG. 3: (Color online) Time-dependent field of the pump (a), 
and dump (b) pulses. The field in each pulse is normalized 
relative to the maximum of the pump field. 

pairs of pump and dump pulses, all shaped as in Fig. [3] 
The trains envelopes vary linearly on the pulse number 
(Fig.jija)). The inter-pair separation is 1310.59 ps, which 
is close to a revival time of the wave packet |2). The car- 
rier frequencies are chosen to match the Raman condition 
w Raman = ujp—ujp,. The carrier optical phase is kept con- 
stant throughout each pulse train by choosing the carrier 
frequency to coincide with a comb tooth [BJ: 

Up.D = ^(Np^ofrep + foP,o) , (2) 

where f rep is the common repetition rate and fo p.o is 
the carrier-envelope offset frequency in the pump (dump) 
train. The process is interferometrically sensitive to the 
inter-pulse delay (f rep ) and requires phase locking of the 
pump and dump combs; however, phase stabilization of 
frequency combs is now a standard technique [5]. 

If the spontaneous decay of the A 1 Y.-b 3 H states is ne- 
glected, the calculation predicts transferring 92% of pop- 
ulation from |1) into |3). Introducing a 15-ns A^-b^Ii 
decay time reduces the transfer efficiency to 75%. The 
maximum combined transient population of all the A l Yi- 
6 3 n levels reaches ^0.04. The transfer efficiency is found 
to be rather insensitive to the field strength profiles in the 
pulse trains, and varies within about 15% if the number 
of pulses is varied between 35 and 300. 

Figure |4jc,d) shows an example of CRP-like PAP. The 
field consists of 20 pairs of pump and dump pulses. 
The pump and dump pulse train envelopes are Gaussian 
(Fig.Qc)). The inter-pair separation is 1309.93 ps. The 
central carrier phase in each field evolves with the pulse 
number according to Eq.([T|), with ap = a.p> = 0.09. The 
calculation with decaying A^-b^Ii levels predicts trans- 
ferring 60% of the initial population into the target state; 
neglecting the decay raises this number to 65%. Approx- 
imately 2% of the population ends up distributed among 
eigenstates other than |3) of the target A 1 S + manifold; 
22% of the initial population are transferred to neigh- 
boring states of |1) in the a 3 S + manifold at the end of 
the process. Unlike in the simple 3-lcvcl case, CRP-like 
transfer in the molecule is found to be significantly more 
sensitive than the STIRAP-likc transfer to the number of 
pulses and the intensity profiles of the pump and dump 
pulse trains. We note that the required train envelope 
modulations (either in amplitude or phase) can be read- 
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FIG. 4: (Color online) Dynamics of the STIRAP-like (a,b) 
and CRP-like (c,d) PAP in KRb. (a,c): Envelopes of the 
pump (dashed blue) and dump (red) pulses. Each point shows 
the factor by which the field shown in Fig. [3] is scaled. The 
carrier phase (f>^ in the piecewise CRP, same for the pump 
and dump, is shown by the gray staircase in (c). (b,d): Popu- 
lation dynamics of all the levels in the system. Time is given 
at the "compressed" set of points: output is shown only for 
the times when the field is on. 



ily generated experimentally by present modulators for 
standard 1-10 ns repetition-period pulse sources. 

The dependence of the transfer efficiency on AT in 
both STIRAP-like and CRP-like PAP is irregular due 
to the complex dynamics of the A^-b^Ii wave packet. 
Viewed from the spectral perspective, each of the A 1 ^- 
b 3 H eigenstates can serve as an intermediate level for a 
separate AP pathway. Reminiscent of the traditional co- 
herent control via a wave packet [U [T5] , the AP pathways 
interfere to provide the population transfer into the tar- 
get state. Thus the frequency comb source drives a co- 
herently controlled AP process [5J [5] via the set of inter- 
mediate levels. Matching the repetition time and carrier 
phase to the wave packet revival time and phase ensures 
that all the possible AP pathways participate; spectral 
shaping of the pump and dump pulses sets constructive 
interference among the pathways. 

Tailoring the optimal pulses requires detailed knowl- 
edge of the excited state structure, which is not always 
available. But what happens if the pulses are not shaped 
at all? Fig. [5ja) investigates the total efficiency of the 
STIRAP-like transfer under both pump and dump trains 
consisting of 50 sin 2 crf-shaped pulses, with respect to 
the inter-pair delay AT — 1/ f rep and intra-pair delay 
ST between the pump and dump pulses. Raman res- 
onance is maintained throughout the f rep scan by vary- 
ing the carrier-envelope frequency difference according to 
fop - foD = (uRaman/Zn) - N f rep . The plot shows a 
remarkable landscape of vertical and horizontal features, 
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FIG. 5: (Color online)(a): An efficiency density plot for the 
STIRAP-like PAP with unshaped laser pulses in dependence 
on the inter-pair delay AT and the intra-pair delay ST. Pos- 
itive ST indicates that the pump comes after the dump; i.e. 
population remains in the excited wave packet until the next 
dump comes (suffering more spontaneous emission losses). 

(b) : Vertical at AT = 1310.62 ps (blue dotted line in (a)). 

(c) : Fourier transform of (b). The green dotted line marks 
the beat frequency between and E 2 . 

where some peaks reach 80% transfer efficiency. Note 
that no population is transferred to any of the a 3 E + or 
AT 1 E + states, except |1) and |3), or to any of the inter- 
mediate states. In the spectral domain, either one state 
guides one of the interfering AP pathways, or it is simply 
missed by the narrow spectral selectivity of the frequency 
combs. In the time domain, destructive Ramsey- type in- 
terference prevents transitions to undesired levels. 

To elucidate some of the structure observed in Fig. [5] 
consider ST expressed in frequency as a linear spec- 
tral phase of the dump comb with respect to the pump 
comb. Accordingly, the phase of a specific Raman path 
at detuning f2j from the carrier frequency is <p(Cli) = 
—iliST + <fid(Qi), where (/>d(f2j) is the relative phase of 
the dipoles associated with this path. The vertical lines 
in the image correspond therefore to a AT value where 
one or several teeth are in one-photon resonance with Ra- 
man paths. Changing ST at this AT reveals interference 
between these paths. A Fourier transform with respect 
to ST reveals the spacing between the intermediate lev- 
els, as shown in Fig. f5Vb,c), where mostly two Raman 
paths, via = 11145 cm" 1 and = 11190 cm" 1 , 
contribute, and the Fourier transform of the transfer ef- 
ficiency peaks at the beat frequency between E% and 

(2) 

E2 ■ Horizontal lines, on the other hand, occur when 
ST hits a special value, where y>d(f2<) ~ SliST for the 
participating Raman transitions; i.e., ST is matched to 
the vibrational dynamics. In this case the transfer peaks 
even if there is no coherent accumulation in the inter- 
mediate levels, as was observed already in [7]. Since the 
match to the dynamics is only approximate, the lines are 
not perfectly horizontal. In the best case scenario every 
dump pulse transfers into |3) exactly all the population 
excited to |2) from |1) by the preceding pump [7]. 



With some preliminary knowledge of the initial and 
target state's energies, the range of frequencies and the 
field strengths needed (but not of the exact values of in- 
termediate energies and transition dipoles), one can take 
trains of unshaped pump and dump pulses, and scan AT 
and ST to find efficient PAP transfers. In an analogy 
with 2D Fourier spectroscopy [T3], a theoretical analysis 
of the 2D scan of the transfer efficiency may be able to 
provide the full spectroscopic information. Further, one 
can experimentally optimize the pulse shapes in order to 
maximize the transfer fidelity. 
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